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Abstract 

A model is proposed which can be regarded as a mean field approximation 
for pure lattice QCD and chiral field. It always possesses a phase transition 
between a strong-coupling phase (where it reduces to a one-plaquette inte- 
gral) and a non-trivial weak-coupling one. For the U(A^) gauge group, it is 
equivalent to some hermitian multi-matrix model. This analogy allows for 
determining possible large critical regimes thus generalizing the Gross- 
Witten phase transition in the one-plaquette model. 



1 Introduction 



The recent few years have seen the considerable development of the planar- 
diagram technique connected mostly with the matrix models of 2D gravity. 
Unfortunately, it seems that all those achievements has brought no new in- 
sights into large gauge theory, which was the original motivation for the 
method Nevertheless, as was demonstrated by the exact solution of 
QCD2 on a sphere 0, a reduction to a hermitian matrix model can be very 
profitable technically. 

In the present paper I establish a connection between a special class 
of hermitian multi-matrix models and a mean field approximation for pure 
lattice QCD. It enables for using the saddle-point technique for the former 
in the rather new framework. The most interesting phenomenon here is, 
probably, a large phase transition of the Gross- Witten type 0, which 
should have some stringy interpretation. For the one-plaquette model, the 
connection with 2D gravity was discussed in a number of papers 0] . 

The standard QCD mean field (MF), as was proposed by K.Wilson p (for 
review see P), suffers from the obvious drawback of being gauge dependent. 
I suggest a purely geometrical approach to the problem, which avoids the 
step of gauge fixing. This model, however, shares all limitations of any MF 
approximation at the price of being, in principle, soluble. One can say it 
learns nothing about QCD itself. Unfortunately, the same could be said 
about any other model solved so far. 

The starting point is to substitute a regular D-dimensional lattice by 
an infinite (Cayley) tree constructed of two-dimensional plaquettes. The 
plaquettes are glued along their edges so that the tree is a simply connected 
covering of the lattice. Therefore, gauge theory defined on such a tree can 
be regarded as a MF approximation for lattice QCD in D dimensions. 

A similar idea was put forward in Ref. ^ where gauge theory on a Cayley 
tree made of cubes was considered. As far as phase structures of lattice mod- 
els are concerned, the cube-made Cayley tree might provide better accuracy 
than the plaquette-made one (although it is not clear a priori). However, 
the latter enjoys the property of being soluble in the large A^ limit by a 
saddle-point technique, while the former can hardly be handled for continu- 
ous gauge groups. Another nice feature of the model under consideration is 
that it includes both chiral field and gauge theory on equal footing , actually 
interpolating between them. 



1 



The transfer-matrix formaUsm is very simple in our case. Let us consider 
a tree with one root, i.e., to leave one of gauge variables, m, not integrated. 
Then, the corresponding partition function Iv{u) for the rooted tree of a 
large volume V obeys the equation 

p p p 

IpmV+i{tJ') ^ l[dxk K{ullxk) n Ivi^k) (1) 

k=l k k=l 

where the tree is assumed to be made of {p + l)-sided polygons, m + 1 on 
each hnk. If p = 3, it is a covering of the D — dimensional hypcr- 
cubic lattice. The value p — 1 corresponds to the case when a plaquette- 
made tree degenerates into an ordinary one constructed of one-dimensional 
links. Hence, we have a model interpolating between the spin and gauge MF 
approximations. For the SU{N) group, it interpolates between Dqf = 
dimensional chiral field and Dqt = Dqf + 1 lattice gauge theory. 

The Boltzmann weight, K{x), is a real positive class function. Two stan- 
dard choices of it are: the Wilson one 

Kix) — exp — -tr ix -\- x"*") (2) 

and the heat-kernel 

K{x)^Y.dre-^''^Xr{x) (3) 

r 

where is a bare gauge coupling; Xr{x) is a character of an irrep C^j. is a 
second Casimir; dr — Xr{I) is the dimension of r. 

In the thermodynamical limit, the conveniently normalized quantity 



J(m) = lim e"(^+^)^/y(M) (4) 
(where / is a free energy per volume) obeys the equation 

J{u) ^Jfl dxk K{u f[ Xk) n J'^ixk) (5) 



k=l k=l k=l 

and the free energy is given by the relation 

2^ _ pm- 1 
m -|- 



m ^ log I duJ-^\u) (6) 
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Eq. (Kp always has the trivial solution 



J. 



3o 



dx K{x) 



mp— 1 



(7) 



In this case, the free energy coincides with the one-plaquette-model one 

fscig') = log / dx K{x) (8) 



This is the strong-coupling phase of the model. 

In the weak-coupling phase, J{x) is a non-trivial class function: 



X) 



(9) 



where the sum runs over all irreps of a gauge group. Eq. (|^) can be rewritten 
in terms of the Fourier coefficients jr as 



jr = K 



1 

dr 



E li js, / dx Xr{x) n Xs,^ 



X] 



(10) 



k=l 



"r si,...,Sm k=l 

where A,, are Fourier coefficients of the Boltzmann weight. 

The free energy can be rewritten as the single sum over representations 



pm-1 V- , . 

- log 2^ drjr 

-L ~, 



m 



Jr_ 



(11) 



A glueball spectrum can be determined from the eigenvalue problem: 



(12) 



where the analog of the transfer matrix in our case is simply 



m 



rt 



J J2 li 3su I dx Xr{x) n XsS^) 



"r si,...,Sm k=l 
m—1 



p-l 



k=l 



d 



m—1 „ m—1 

r si,...,Sm-i k=l k=l 



(13) 



X) 



Eq. (12) always has the trivial eigenvalue rric = 0, which corresponds to 
the identity operator {i.e., the partition function). The number of possible 
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excitations in the system is equal to the number of irreducible representations 
of a gauge group. In the weak-coupling phase they are all exited. However, 
non of them can become massless and, hence, there is no continuum limit 
associated with the model. 



2 Phase transition 

Let us assume that, at some critical value g"^, the Fourier coefficients of all 
non-trivial representations in Eq. (^) vanish {i.e., the weak-coupling solution 
transforms smoothly into the strong-coupling one). Then, in the vicinity, 
Ag"^ = g1 — g"^ <^ 1, the most essential contribution comes from the funda- 
mental representation, r = 1, 

Jo = 0(1); n = 0{/\g'')- jr = o{Ag'), for r ^ 0, 1 (14) 
After expanding Eq. (p!OD in ji, one finds 

n = xi[2jr'jir+--- ^^^^ 

From which it follows that, for p 1, the smooth transition between the 
strong and weak-coupling solutions is impossible. For p = 1, the critical 
value gl is determined by the equation 

1 ^ AML (16) 

m diXo{g'i) 

For SU{N) with the Wilson weight, one finds 

1 = l(tr/7) = \ (17) 

Hence, the phase transition takes place when the one-plaquette model is in 
its strong-coupling regime [i.e., above the Gross- Witten critical point). Also 
notice that, in the m oo limit, the strong-coupling phase disappears. For 
the heat-kernel, one finds simply g"^ = logm independent of A^. 



4 



The lowest two representations play a special role in the p = 1 case. 
Therefore, the simplest Z2 model may be rather instructive as far as the 
nature of the phase transition is concerned. In this case one can represent 

J{x)=Jo+JiX = pe^'l' (18) 
where x = ±1 is a Z2 variable and 

jo = pcoshcf) ji = psinhcj) (19) 
Let me choose Aq = 1, then Eq. (p!0D takes the form 

p cosh (/) = {p^ cosh uKpY 

p sinh (j) = Xi (p™ sinh riKpY 



(20) 



or, equivalent ly. 



, , , .VI ( cosh \ mp-1 

tanh = Ai tanh'^ m0 p = I — — - I (21) 



Vcosh^ nKp^ 

Hence, a solution can always be determined from some algebraic equation 
with respect to x = tanh (j). 

For example, for m = 3 and p = 1 (2D Ising), one finds the equation 

, 3x + 

from which 



3Ai - 1 1 8Af 

" = ^/^^ ^=2^"g 6A,-Af-l 



And the transition is of the second order: /'(|) = 0. This is the case, if 
p = 1, for arbitrary m and all compact groups. In general, it can be easily 
proven by using the representation (|Tl|) for the free energy. 

In the Z2 model, there is only one excitation for which, in the previous 
example (m = 3, p = 1), one finds the mass nic = log (jz^p- One can easily 
show that, at the critical point for all m, mG(A^) = logm. Therefore, the 
transition does not produce any continuum limit. 

Of course, when p = 1, one just repeats the standard spin MF approxima- 
tion H]. If p 7^ 1, the phase transition between the strong and weak-coupling 
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phases is of the first order. Eq. ( pi]) has obviously no real solutions for Ai 
small enough. Hence, the weak coupling branch disappears somewhere being 
already meta-stable. The first-order transition point can be determined from 
the equation 

p+i _i p+i 
Jo" =1 (24) 

which is quite a standard numerical problem. 



3 Connection with mult i- matrix models 

The Fourier representation (|TU]) is not convenient for the investigation of the 
weak-coupling phase. In this case, the original matrix variables are more 
suitable. Let us choose the Boltzmann weight in the form of the U{N) heat- 
kernel: 

K{x) = Y.dre-^''^Xr{x) (25) 

r 

where, in terms of the highest weight components (mi > m2 > ■ ■ ■ > rriff), 

Cr = Y.{^k-k + N) (26) 

k=l 

is a conveniently redefined second Casimir. For diagonal matrices, x = e*", 
one has the Weyl formula for characters 

- ^ (27) 

where 



AJe'") = det e*('»fe-fc+^)". 
A(e^") = det e^^^-'^)"^ = Y[{e'''' - e'"^) 



(28) 



]<k 



^. = X.(/) = n(l + ^^^^^) (29) 



k- j 
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is the dimension of an irrep r. 

Let us substitute (|^) in (H) and integrate over angular parts of all vari- 
ables: 

1=1 k=l r k=l 



i=l 
N-l 



(n-0 



p-i 1 



2 N 



n=l 



n 



2 N N 



N-l 

( n "0" 



fc=l 



n E 



, nfe=— CO 



dAfc [A(A)] 



exp 



y- V- X / xl /r dete^^^-^»* 

fc=l fe=l 1=1 ^ 



The last equality holds owing to Poisson's formula. Now, Eq. (^) takes the 
form (after renormalizing J{x)) 



AT 



k=l 



N / +00 

n E 

fc=l \nfe=— oo 



dXk A (A) exp 



k=l 







HP 



(30) 



fe=i 



A(A) 



Let us introduce a new function -F(a) such that 



J2 A(a + 27m)F{a + 27m) 



(31) 



which obeys the equation 



F{a) 



A(a) 



N 



k=l 



2aT N 



N 



nrfAfcA(A)exp -^^A^ + ziV^Afcafc 



k=l 



k=l 



+ 00 



N 



k=l 



J™-i(e*^)F(/3) 



(32) 
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If one introduces hermitian matrices A, A, B having eigenvalues a, A, /5 
correspondingly, then one can rewrite Eq. (^) as the matrix integral equation 



F{A) 



^Af2^ g-ig2;vtr A^+iATtrAA 



'.Ntl AB+V\B] 



F{B) 



= (-) ff[d^'B.e-^'''^^^'^''^^'^^'^^f[FiB.) 
\ 9 J i=i i=i 

where the effective potential is determined by the equation 

y[S] = (m-l)log E ^^^^^W + 2vrn) 



(33) 



(34) 



V[B] is a symmetric non-singular function of eigenvalues, hence, can be ex- 
panded in Schur polynomials. It is easily checked that, if m = 1, F{(3) is 
Gaussian and J(e*^) is given by the heat-kernel in the Dowker form |p. 



4 Large N solution 



If p = 1, nothing prevents one from introducing an arbitrary potential, 
NtTU{B), for the chiral field. In this case, what one ends up with is a 
unitary analog of the Bethe-tree matrix model considered in Ref. [|10|- One 



can use a similar technique in both cases. The first step is to write down the 



self- consistency equation for the resolvent matrix [11 : 



F{A) 




-tr {A+B)^+V[B]-NtT U(B) 



F{B) 




N 



A 9^ d 



1 



NdA^ z-B 



^~^il{A+Bf+V[B]-NtlU{B)p^^^ 



N t^-. \ N dak N ak- aj 



k=l 



(35) 
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where 



G,{z) 



'2n 
^ 9 





1 


1 dB 




_z-B 



^-J^^tT {A+B)^+VlB]-NtTU{B) ^^^^ 



(36) 



J kk 



are diagonal matrix elements of the resolvent matrix. This equation is a 
recursive relation for moments of B. 
Let us introduce two functions 



F(x,z) 



, 1 

TT^r — 

^A^ X 



— tr 

X- Az- 



dy 



p{y) 



X 



y 



-X 



dy 



p{y)W{y,z) 



(37) 



B J ' X — y 

p{y) is a density of eigenvalues; W{y, z) is a real function on a support of 
p{y). As we are looking for a homogeneous ground state, /(x) is the same at 
all sites of the tree and F{x, z) is symmetric: F{x, z) = F{z, x). 

From Eq. (^) it follows that, at = oo, W{x, z) obeys the equation 



{z + g'w{x) + x)W{x, z) + g' I dy p(j/) ^^'''^^ W{y,z) ^ ^ 

X — y 



where 



1 d 

w(x) = lim — — — logF[y41 



d 6 



lim 



1 



0'k=x dx 6p{x) N^oo A^2 



logF 



(39) 



As was first noticed in Ref. |T^, equations of this type can be solved by 
the Riemann-Hilbert method. The outcome of which is the following integral 
representation 



F{x, z) = 1 — exp J 
where, in our case. 



dy 1 z-u+{y) 
log- 



27ri X — y z ~ U-(y) 



(40) 



u±{x) = —X — g'^{w{x) + Re /(x) ± z Im f{x)) 



(41) 
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and the integral goes along a support of Im f{y) = np{y). Off the support of 
p{x), Ref{x) and Im/(x) continue analytically as two independent holomor- 
phic functions. Eq. (^) makes sense as a set of recursive integral relations 



obtained by expanding both sides in inverse powers of x and z (see Ref. [|T0 
for details). 

As F{x, z) is symmetric, the following equation holds 

u+{u-{x)) = X (42) 
The function w{x) can be determined from the saddle-point equation 

2Re/(x) + 2w{x) + Jhn ±V[x] - U\x) = (43) 

which is, in general, a complicated non-linear integral equation. However, for 
small enough, we can neglect in the large limit all non-trivial winding 
numbers in Eq. (0) {i.e., put Uk = \/k), then Eq. (|^) takes the simple 
form 



2Re f{x) + {m + l)w{x) — {m. — l) J dy p{y)(- cot 



x-y 




or 



1 f 1 1 

w{x) = ^-\u'{x)-2Ref{x) + {m-l) Y.'[f{x-27tn) + —-]\ (45) 

+ I n=-oo J 

The prime means that the n = term in the sum is omitted. This formula 
makes sense only when a support of p{x) lies inside the interval (— tt, +7r), 
i.e., when there is a gap in the eigenvalue distribution of original unitary 
matrices. 

If p > 1, any local potential for unitary variables spoils the gauge invari- 
ance of the model. However, one can introduce a non-gaussian potential for 
the auxiliary field A in Eq. (|33|). It corresponds to taking an arbitrary Boltz- 
mann weight. In this case, one cannot easily integrate out the auxiliary field. 
Nevertheless, our method can be generalized to the inhomogeneous system. 
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Now, one has two different external field problems and has to introduce three 
functions 



viz) = — tr = 



z-X 



^/ X , 1 1 1 ^ 

F(x,z) = — tr -) = 



dy 



p{y)W{y,z) 
x-y 



dX 



r]{X)n{X,x) 
Z-X 



(46) 

Both W{x, z) and ^(A, x) obey equations analogous to ( ^Sf ) and -F(x, z) has 
two different integral representations of the type (0). For example, 



I{X)=JdB e*^tr^^+^[^]F(5) 

dB—ti (xH — 

J N V NdA^x-1 



,iNtT AB+V[B] 



F{B) 



(47) 



TV 



i d 



fc=i 



G',( a:)-G',(x) 



where Gk{x) are diagonal elements of the resolvent matrix as in Eq. 
Then one finds 

{x + zuj{X)MX, x) + i [dfi ^(^) ^(^'^)-^(^^^) ^ 1 



where 



X — 

d 6 



(48) 



Afe=A aA 5r7(A) ^-oo A^2 



hm —/(A) (49) 



The saddle-point equation with respect to A gives 

2Re ifiX) - U'{X) + {p + l)u{X) = 
where U{X) is an arbitrary, in principle, potential. 



(50) 
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Thus, one finds for F{x, z) tlie representation 

F{x, z) = l-exp f ^ log - — ^^iy^ 



wfiere 



v±{\) = -{^U'iX) + ^Re^(A) ± tlm^iX)} 
In a close analogy, one finds 

{z + iw{x))W{x, z) + i[ dy p(^y^ ^(^^^)-^(y^^) ^ 1 

J X — y 

and the saddle-point equation for a distribution with a gap: 



(51) 



(52) 



(53) 



2Re/(x) + (m + l)w{x) - (m - 1) J2'[fi^ " 2™) + 2^ 

n=— 00 

Hence, 

F{x, = 1 - exp / — log -- 

J Imx — y z — U-{y) 

where 



(54) 



(55) 



«±(x) = -{^-Re/(x) + ^- J2'[f{x~2nn) + ^]±zlmf{x)} (56) 



And, instead of Eq. (^2]), there are two equations 

ti+(f_(x)) = X w_(t;+(a;)) = x 



(57) 



5 Critical regimes 



Following Ref. one can determine possible large critical regimes. 

Eqs. (1^) and (^) allows, in principle, for constructing the functions f{x) 
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and if{x) thus solving the model at large A^. Practically, it is a very com- 
plicated problem. However, universal behavior is determined by a scaling of 
the imaginary parts of f{x) and (p{x) near their edges, which are, in general, 
branching points: 



f{x) = freg{x) + c{x - Xq) 



1+7 



^{\) = ^reg{\)+c{\-\of^^ ... (58) 



/reg(x) and (preg^x) are regular parts of the functions at the branching points. 

Let us start with the p = 1 model. We are interested in the situation 
when the edges of the distribution of eigenvalues of original unitary matrices 
collide. It corresponds to the case when the edges of Im/(a;) in Eq. (^) 
touch branching points of /(a;±27r). Let us expand all quantities in Eq. (p]) 
near one of the collision points (by redefining the variables one can always 
place it at the origin): 



u±{x) = ax + b 



m 



m + 1 



(cos7r7 + l)e-'''^x^+^ ± ib sin nje-''''^x^+^ + . . 



(59) 



By substituting it in Eq. (^21) , one obtains two equations 



and 



m 



m+ 1 



[cos 7r7 + 1) — 2 sin + a"' 



ym — 1 
m + 1 



(60) 



cos 7r7 + 1) + i sin 7r7 =0 (61) 



The simplest possibility is a = 1, when one finds 



cos 7r7 



-1 



(62) 



which is the standard Gross- Witten singularity always possible in the model. 
By tuning the potential U{x), one can reach a multi-critical point when 
a = — 1, then there are two branches: (i) a* = e^'^^ yielding the equation for 

7 

1 

cos 7r7 = — 
ni 
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(63) 



and (ii) a? = e "^'^ yielding 

cos 7r7 = m (64) 

Let me remind a reader that m is connected with the effective space 
dimension as m = 2D — 1. The obvious duahty m — is iust a manifestation 
of the duahty D 2D-1 Bethe-tree matrix model found in Ref. [|1^ . 

The solution fits two soluble cases: m = — 1, which corresponds to the one- 
plaquette model, and m=l, which is just a one-dimensional matrix chain. 
The case m = corresponds to a two-matrix model (as we have introduced 
an arbitrary potential, it is not simply reducible to a one-matrix integral any 
more) . 

li p > 1, there are two functions obeying Eq. (0). Expanding them as 

777- — 1 

u±ix) = ax - b (cos7r7 + l)e"'"''^a;^+'^ ± ib sin n'ye~'''^^x^~^^ + . . . 

m + 1 

v±(x) = a'x - cos7r7e-*"^x^+^ ± sin ^7e-*"^x^+^ + . . . 

^ ^ p+l 

(65) 

and substituting the expansions in Eqs. (0), one finds that aa' = 1 and 
a'b (cos7r7 + l)±isin7r7 —b'a^'^'^ cos 7r7 ± i sin 7r7 =0 (66) 



m + 1 



p+l 



This is a degenerate system of linear equations for the real and imaginary 
parts. Equating to zero the corresponding determinant, one finds 

cos7r7 = - (67) 

~'~ p+l m-l 

At p = 1, this formula reproduces Eq. ( p^ ) and, at m = 1, gives 7 = |, which 
is the simplest allowed singularity in a general matrix model. 

The critical regime ( p7D corresponds to the situation when the potential 
for the A variable becomes critical (in the standard matrix-model sense) at 
the same moment as the edges of the density for the unitary matrices collide. 
It is a multi-critical point with respect to the Gross- Witten singularity, which 
is always allowed and corresponds to the case where M+(a;) has a quadratic 
extremum matching with a square- root branching point of z;_(x). 
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6 Discussion 



Our method allows for determining possible large critical regimes without 
really solving a model. It means that some critical points can correspond to 
non-stable or non-unitary models. It was indeed the case for the Bethe-tree 
matrix model at D > 1, where critical potentials appeared to be, in general, 



complex |T0|. As our MF approximation for chiral field is quite similar to 
the one for hermitian field, one can expect that the m > 1 branch (|6^) 
for the p = 1 model exists only owing to the formal duality, in the sense 
of analytical continuation into an unphysical region of parameters. On the 
other hand, the scaling (^) in the gauge MF model should be quite sensible 
for p > 1 and m > 1 as corresponds to the simplest reachable singularity. 
However, to determine a critical form of the Boltzmann weight is the crux of 
our approach. 

Usually, the Gross- Witten transition is considered as a pure lattice arti- 
fact having no physical meaning. However, it is a very general phenomenon 
taking place in any model reducible to a saddle-point problem for a unitary- 
matrix-valued master field. 
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